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ABSTRACT: 


We present here the results of detailed Monte Carlo calculations of the 
interaction histories of ultrahigh energy cosmic-ray nuclei with intergalactic 
radiation fields using improved estimates of these fields and empirical determ- 
inations of photonuclear cross sections including multinuclear disintegrations 
for nuclei up to ^Fe . Intergalactic and galactic energy loss rates and nucleon 
loss rates for nuclei up to Fe are also given. Astrophysical implications are 

discussed in terms of expected features in the cosmic-ray spectrum between 10 18 
21 

and 10 eV for the universal and supercluster origin hypotheses. The results of 
these calculations indicate that ultrahigh energy cosmic rays cannot be universal 
in origin regardless of whether they are protons or nuclei . Both the supercluster 
and galactic origin hypotheses, however, are possible regardless of nuclear com- 
position. 
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1 . Introduction . 

The basic interactions between photons and nuclei of ultrahigh energy which 
have important astrophysical consequences are: 

1) Compton interactions, 

2) Pair production (particulady of e + e pairs) in the field of the nucleus, 

3) Photodisintegration of the nucleus, 

4) photoproduction of hadrons 

in the rest-system of the nucleus, process (1) has no threshhold energy, 

2 

process (2) occurs at a threshhold energy of 2m e c ~ 1 MeV, process (3) is 
particularly important in most cases between 15 and 25 MeV where the giant dipole 
resonance has its peak, and process (4) has a threshhold energy of 

e ' - m c 2 (l+m /2m_) - 145 MeV (1) 

th TT TT P 

Process (1) , Compton effect interactions involving nuclei, is a process 
which results in only a negligibly small energy loss for the nuclei given by 



-3 

where P-, is the energy density of the ambient photon field in eV cm , a T 
is the Thomson cross section and Z and A are the atonic number and weight of 
the nucleus with E being the total energy of the nucleus in eV. 
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Processes (2) through (4) must be considered in investigating the propa- 
gation ultrahigh energy cosmic ray nuclei in intergalactic space (Greisen 1966, 
Zatsepin and Kuz'min 1966, Stecker 1968,1969, Berezinskii and Zatsepin 1971)’. 

Process (2) , particularly important for heavy nuclei, has been treated in 
detail in its astrophysical context by Blumenthal (1970) . This process is only 
of importance for interactions with the 2.7 K blackbody background radiation 
(see section 4). In this case, for interactions with a blackbody field of temp- 
erature T, the photon density is given by the Planck function 

n(e) = (Hc)- 3 (eA) 2 ( e E/kT - 1 j" 1 (3) 

and the energy loss rate for the nucleus is 

-(dE/dt) = 3ao m h -3 Z 2 (m c 2 kT) 2 f {m c 2 / (2ykT) } (4 * 

T s 6 

/ 2,-1 

where f ( v ) is shown in Figure l.The Lorentz factor of the nucleus y = EfAm^c ) . 

The photohadron production process (4) is dominated by photopion production 
channels. This process has been treated in detail in a previous paper (Stecker 
1968) . Figure 2 shows the lifetime and mean free path for protons against atten- 
uation by energy loss due to pair production (Blumenthal 1970) and photomeson 
production (Stecker 1968) , due to interactions with the 2.7 IC radiation. Also 
shown in Figure 2 is the lifetime for Fe against energy loss by e e pair 
production. 

Process (3) , photodisintegration, is not as easy to treat as the other pro- 
cesses and most previous cosmic ray papers have used only simple estimates of 

4 

the cross section for this process, stecker (1969) used a detailed study of He 

4 

photodisintegration (Gorbunov 1968) to determine the lifetime of He and estimate 
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the lifetime of Fe against one nucleon loss by interactions with the 2.7 K 
radiation. The main emphasis in this paper will be on a presentation of the 
results of a new and much more detailed calculation of the intergalactic 
photodisintegration of ultrahigh energy nuclei with the following improve- 
ments : 

1) Use of empirically determined cross section data as a function of 

energy for all nuclei with l^A-56. It is exactly in this mass range where 

contributions from nonresonant photonuclear disintegration processes are most 

56 

significant. Because of the position of Fe on the binding energy curve, it 
is considered to be a significant end product of stellar evolution and highter 
mass nuclei are found to be much rarer in the cosmic radiation. 

2) Use of cross section data for 15 a£ 56 for multinucleon emission inter- 

actions as well as single nucleon emission interactions. Multinucleon losses 
(except for double nucleon loss) involve non-resonant processes and are particu- 
larly important at energies between the position of the giant dipole resonance 
and the pion production threshhold ( i.e. , 25 MeV<s ' <145 MeV , see section 

2 ) . 

3) Use of new estimates of the intergalactic infrared radiation field, whose 
effect is included along with that of the blackbody radiation. 

4) Inclusion of energy loss by pair production occuring concurrently with 
energy loss by photodisintegration. 

5) Treatment of these Concurrent energy loss processes by using numerical 
Monte Carlo techniques which also allow a natural treatment of the multinucleon 
loss processes in a determination of the photodisintegration and energy loss 



6 


histories of nuclei as a function of initial energy. 

2. The Photodisintegration Process . 

A general discussion of the photonuclear interaction process for the 

energies of interest here has been given by Danos and Fuller (1965) and 
Hayward (1970) . This energy band, extending up to 150 MeV, or about the 
meson production threshold, rather naturally splits into two parts. The lower 
region, extending up to 30 MeV, has been studied extensively. This is the 
domain of the giant resonances of the nuclear photoeffect. Here, although two- 


nucleon emission does take place, most of the absorption cross section results 
in the emission of only single nucleons, either neutrons or protons. For the 
medium and heavy nuclei, A > 50, the total photon absorption cross section can 
be well represented by a single, or in the case of the deformed nuclei, by 
the superposition of two Lorentzian curves of the form 


o (e 1 ) = 


(e'g - 


Td 
1 2 \ 2 


)* + e ,2 Id 


(5) 


For the lighter nuclei, the absorption cross section often shows considerably 
more structure than this. For all nuclei the basic photonuclear interaction at 
these energies is pictured as being with the individual uncdrrelated nucleons 
com p rising the nucleus. In the shell model language, the interaction results 
in the creation of particle-hole pairs. The collective giant resonance state 
is described by a coherent superposition of individual single-particle transi- 
tions . 

The data for energies between 30 and 150 MeV are relatively sparse . The 
available measurements seem to indicate that the total cross section in this 
region is a relatively smooth function of energy. Reactions in which a number 
of nucleons are emitted seem to make up most of the total absorption cross section. 
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In this region the photon's wavelength becomes comparable to, and less than, 
nuclear dimensions and it becomes very difficult to conserve both energy and 
momentum in a single-particle interaction within the nucleus. The dominant 
absorption process is then pictured as one in which the photon interacts with 
a nucleon pair while they are scattering within the nucleus. This is referred to 
as the quasi-deuteron effect. (Levinger, 1951). 

In the following discussion, the strength of a photonuclear reaction in the 
two energy regions, i.e., the integral of the cross section over the energy 
interval, will be given in units of the classical dipole sum rule, i.e.. 



= ja (e ' ) de ' = 


2tr^e^ 

Md 


NZ _ 
A 


59.8 MeV-mb 


( 6 ) 


where M is the nucleon mass, A the mass number, Z the number of protons, and N 
the number of neutrons in the nucleus. The derivation of this result is based 
on both the long wavelength approximation for the dipole operator involving 
nucleon coordinates only and on closure, i.e., it assumes that the integration 
in Eq. (6) is. carried out over all excitation energies. Meson coordinates within 
nuclei are specifically ignored. Since these assumptions are not valid for the 
excitations being considered here, no real physical significance should be attached 
to this normalization. It is simply a convenient scale factor that removes the 
principal dependence of the total absorption cross section on N, Z, and A. 

The data used in these calculations were based on an evaluation of the information 
in the files of the Photonuclear Data Center at the National Bureau of Standards 


(E. G:. Fuller, private communication). To simplify these calculations, the observed 
detailed dependence of the various reaction cross sections on photon energy was 
not used. Rather, each cross section was assumed to be independent of energy over 
a finite energy band of width A . For energies below 30 MeV, this band is 
assumed to be centered about a mean energy e '' for each reaction. As a further 
simplification, it was assumed that the half-lives of the various reaction 
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products resulting from photonuclear reactions were all short compared to the 
rate at which the reactions take place. With this assumption, it was not necessary 
to consider each individual reaction but rather only the cross section for 
emitting AA nucleons from a nucleus of mass A had to be input into the calculations. 

The data used in the calculations are given in Tables 1 and 2. For energies 
below 30 MeV, the total absorption cross section is assumed to result in reactions 
in which only one or two nucleons are emitted. These are given in terms of the 
quantities e Q 1 , £, and A listed in Table 1 with the subscripts 1 or 2, respectively. 
For each reaction, these were defined as follows: 



,30 ( 1 ) 

J a 1 ' (e')de* 
0 


(7) where a^(e') is the actual single nucleon 
emission cross section as a function of photon 
energy e ' . 

O' 



the mean energy for the reaction a^{£' ) chosen by inspection 
of a graph of the cross section as a function of e'. 


,30 /i\ 

A, a = I a (e'Jde' (8) A, and o were chosen so as to closely as 

1 m 4 0 l m 

possible match the actual dependence of 
a^-(e') on energy. 

The corresponding quantities for the two nucleon reactions were chosen in a 
similar way. As has been indicated previously , there are ample experimental 
data available for photon energies below 30 MeV to use in determining these 
parameters. For these excitation energies, the principal undertainties in 
the final results are probably those associated with the very crude parameteri- 
zation used to describe the one and two nucleon emission cross sections as a function 


of energy 
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The data for the region from 30 to 150 MeV are parameterized by two 
quantities, the total strength of the interaction given by 5 in Table 1, 
and a branching ratio given in Table 2. The total, strength is given by: 


150 

5 - -i- f o fc (e ' >de ' 
Ed 3° 


(9) 


where a t (e') is the total photonuclear absorption cross section. The 
strength between 30 and 150 MeV was actually obtained by taking the 
difference between the estimated total strength to 150 MeV and a value for 
the measured strength in the region 0 to 30 MeV. The estimated total 
strength to 150 MeV as a function of Z was based to a large extent on an 
interpolation and extrapolation of the total cross section measurements 
made at Mainz (J. Ahrens et al. , private communication, to be published in 
Nuclear Physics) . These values are estimated to have an overall uncertainty 
of ±15%. The branching ratio data given in Table 2 for the nuclei heavier 
than Na are based on spallation product yields induced by bremsstrahlung 
with peak energies ranging from 150 to 240 MeV while the data for the 
lighter nuclei are based on the analysis of 170 MeV bremsstrahlung yields 
obtained from cloud chamber measurements. These data should not be taken 
to represent anything more than a first approximation of the yields that 
might be expected from a 1/e' photon energy spectrum. 

For the lighter elements it is found that only about half of a "dipole 
sum” as given by equation (6) is in the giant resonance and even for 39 K 
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only about 3/4 of the value given (6) is in the giant resonance. In the 
particular cases of elements with A 1 4 # the giant resonance does not play 
a role in the photodisintegration process (see later discussion) . 

One may also note that in terms of total integrated cross sections 

30 MeV 150 MeV 

/ a (e-)de'- ^ / o (E')de' (10) 

0 s 30 *- 

which is another indication of the effect of nonresonant processes since 
higher order multipole resonances are considerably weaker than the dipole 


resonance 
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Table 1. Cross Section Parameters for Photonuclear Emission Processes (e and A in MeV) 
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Table 1, Continued 
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Table 2. Branching Ratios for i-nucleon Emission in the 30 MeV to 
150 MeV Energy Range. 


Nucleus 

i 

f . 

i 

Main Channel 

4 He 

1 

0.8 

4 He (y,n) 3 He, 4 He ( Y ,p) 3 H 


2 

0.2 

4 He (y,np) 2 H 

9 Be 

1 

1 

9 Be (y,n) 2 4 He 


10„ 22 
B - Ne 

1 

2 

3 

4 

5 

6 

0.1 

0.3 

0.1 

0.1 

0.2 

0.2 

23 56 



Na - Fe 

1 

0.10 


2 

0.35 


3 

0.10 


4 

0.05 


5 

0.15 


6 

0.045 


7 

0.04 


8 

0.035 


9 

0.03 


10 

0.025 


11 

0.02 


12 

0.018 t.-. . 


13 

0.015 


14 

0.013 





3. The Intergalactic Background Radiation Fields. 


The background radiation of importance in determining the photonulear 
reaction rate in intergalactic space consists of the 2.7 K microwave radiation 
.and the radiation fields of higher energy , viz. , the infrared and the optical 
background radiation. The long wavelength radio background from extragalactic 
sources ( Clark, et al. 1970 ) is of too low an energy and intensity to be of 
importance in the context of this paper. 

a. The Microwave Radiation. 

The 2.7 K radiation appears, for all practical purposes here, to be an 
ideal blackbody in the energy range 2 x lo “ 6 eV<e <2 x 10” 3 eV ( Peebles 1971, 
Thaddeus 1972, Hegyi et al. 1974, Robson et al, 1974, Woody et al. 1975 ). We 
assume it to be such in our calculations. 

b. The Optical Radiation . 

In considering the optical and infrared radiation in intergalactic 
space, one must keep in mind that even for measurements made outside the atmo- 
sphere, the measured background contains components from such nonextragalactic 
sources as Zodiacal light, emission from interplanetary and interstellar dust, 
and radiation from nearby stars at high galactic latitudes as well as halo 
stars. Because, therefore, only upper limits are available from direct meas- 
urements, the intergalactic infrared and optical radiation must be estimated 
using theoretical arguments. 
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Measurements of the fluctuations of the observed optical background 
on at small angular scale give ijpper limits on the contribution to the extra- 
galactic flux from discrete sources such as galaxies { Partridge 1974, Shectman 
1974 ) . The implied flux gives a value for the energy density of the optical 

“3 • *3 

radiation of ^ 3 x 10 eV cm J in the energy range 0.8 eV <e < 10 eV. This 
value precludes strong evolution of the luminosity function for galaxies ( i.e., 
galaxies much brighter in the past ) as has been discussed by Partridge and 
Peebles 1967a, b. Of course, a diffuse field of truly extragalactic origin 
( e.g., bremsstrahlung from a hot intergalactic gas ) would not be detectable 
using fluctuation measurements, but we can assume that such radiation is relatively 
small ( Field 1972 ) . 

We therefore assume an intergalactic dilute starlight field with a temp- 

-15 . 

erature of 5000 K and a dilution factor of 1.2 x 10 in our calculations. The 
effect of such a field is small, as previously noted in the case of a somewhat 
larger optical radiation field by Stecker (1969) . 

-3 

c. Infrared Radiation { 2x10 eV to 0.8 eV ) . 

Only experimental upper limits exist in the infrared energy range 
for the extragalactic background radiation { Pipher et al. 1971, Soifer et 
al. 1971, Houck et al. 1972, Meullner and Weiss 1972, Williamson et al 1973, McNutt 
and Feldman 1970, Hof fmann et al. 1973 ). Thus,: we must rely on theoretical esti- 
mates of this radiation. Here, as in the case of the optical background rad- 
iation, we shall assume that the infrared background originates mainly in 
extragalactic objects. In the last four years it has been found that many of 
the active extragalactic objects such as quasars, N galaxies, the nuclei 
of Seyfert galaxies, and some radio galaxies as well, emit most of their energy 
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in the far infrared in the 30 to 100 u range. Neugebauer et al (1971) , Rieke and 

/ 

Low (1972) and Kleinman and Wright (1974) give fluxes at 10y for 64 of these 
sources. For two sources, M82, and NGC 253, a continuous spctrum has been meas- 
ured. Using these data, combined with the data of Schmidt (1971) on the luminos- 
ity function of extragalactic objects ; it is possible to evaluate the contribution 
of these extragalactic objects to the infrared background. The data indicate that 
a typical extragalactic infrared source has a maximum in its spectrum between 30y 
and lOOy and that the average spectral index of its energy spectrum on the high 
frequency side of the maximum is a - 1.5. 

The effect of sources at high redshifts z can be estimated as follows: 

If one considers class of objects with a density evolution with redshift propor- 
tional to ( 1 + z ) m ( m = 3 being the case with no evolution and only geometrical 
expansion included )and we consider the part of their spectrum which has a power 
law form so that their intensity is given by 

I(v) = v'“ for v <v < » (11) 

0 m M 

their integrated contribution to the infrared background is then 

E (v) = n 0 I 0 v'“ (cH5 1 /4tt) | M dz (l+z) m_ (a+5> (l+ttz)“ 1/2 (12) 

where H Q is the Hubble constant and Q is the fraction of the critical density. 

The solution to equation (12) will give power law spectra in the 
frequency ranges (l+Zj^T^ v < v m and v m £ v £ v M with a spectral index 
depending on the values of a and m. 

If most quasars have a high infrared output compared with their output 
at other wavelengths, it then seems to be a likely possibility that there is a 
strong redshift evolution of their infrared luminosity as is found to be the 
case at optical wavelengths ( m = 6 ) . Considering the limited amount of observa- 
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tional data on Seyfert galaxies, both strong evolution (m=6) and no evolution 
(m=3) have been considered here in estimating their contribution to the infra- 
red background. The fluxes from quasars and Seyfert galaxies, as well as the 
flux from normal galaxies with m=3, have been computed and the results of these 
calculations are shown in Figure 3. 

Two power laws have been drawn on Figure 3 which can be considered as 
giving reasonable lower and upper limits on the intergalactic infrared flux. 

At very long wavelengths (^500vi) the major contribution to the infrared back- 
ground is probably due to the diffuse component associated with the distortion 
of the 2.7 K blackbody spectrum due to effects associated with the early thermal 
history of the universe ( KOmpaneets 1957, Zel'dovich et a! 1972, Chan and 
Jones 1975 ) . 

The power law spectra shown in Figure 3, which were employed in this 
paper , were taken to extend from the optical graybody radiation field to the 
microwave blackbody radiation field , will hereafter be designated as HIR 
{high infrared case) and LIR (low infrared case. They are given in terms of 
photon density spectra by the numerical expressions 

HIR: n(E) = 1.1 x 10~ 3 e 2-5 (cm” 3 eV _1 ) 

and 

LIR: n(e) = 2.6 x 10 -3 e” 2 (cnf 3 eV 1 ) 

with the photon energies given in eV and the spectra extending over the energy 
rancje between 2 x 10~ 3 eV and 0.8 eV. 
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4« Calculation of the Intergalactic Nucleon and Energy Loss Rates . 

The photodisintegration rate for a nucleus of mass A with the sub- 
sequent release of i nucleons is given by the expression 

00 2y_ e 

R. . = — ■ y 2 fde e 2 n(e) / de' e' a . (e') (14) 

A,i 2 'A q J Q A,i 

(Stecker 1969) where n(e) is the photon density of the ambient radiation 
with energy e in the observers system (see previous section), e" is the 
energy of the photon in the rest system of the cosmic ray nucleus and y^ is the 
Lorentz factor of the nucleus. 

Because, as was discussed in section 2, there can exist large and 
significant deviations from the pure Lorentzian shape in the energy range 
where the maximum cross section occurs, a Gaussian approximation for the 
cross section for 1 and 2 nucleon disintegrations in this energy range was 
found to be useful both as an adequate fit to the cross section data and as 
an expedient in performing the numerical calculations. This choice, in eon- 
juction with a threshhold energy of 2 MeV in the rest system of the nucleus 
and a constant cross section above 30 MeV, was numerically found to give 
rates which were very similar to those given by a Lorentian over the resonance 
region where integrated over the photon density function used in equation (14) 
as discussed in the previous section. The error introduced by this approximation 
was found to be relatively small compared to the other uncertainties involved 
in this calculation. Relatively more important, was the consideration of multi- 
nucleon emission channels which, when weighted by the number is nucleons emitted, 
can be seen to have an important effect of the total nucleon and energy loss 
rates. This can be seen from an examination of the branching ratios for multi- 
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nucleon emission given in table 2. 


The cross sections used for evaluating equation (14) were of the form 


* 


K ‘“'Vi 5Ji ^ 




(IS) 


for i m 1, 2 , and 


o- . sftigm 
*■' 110 


for i > 2 


(16) 


where all the energies are in MeV and the normalization factor 


with 


and 


w 


■(!)'[« (&')“'<&)] 


Erf (x) = Ztr 






I 

o e'<x 


6 (*) 


( o e 

T, c' 


< x 


(17) 


(18) 


(19) 


By substituting equations (15) through (19) in equation (14) , the 
photodisintegration rates can be expressed as sums of integrals of four 
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first is 


r 

( &( t *-')■' Jfol 


h/2/ J 


for integrations of the Gaussian part of the cross section over a thermal 

<• ' ‘ • . 

photon field where 


C, = ( ir 2 R 3 c 3 ^CVT 1 ? .ZJT 1 
A A,i d 


D being a dilution factor for the case of starlight optical radiation (graybody) 
which is equal to unity for pure blackbody radiation* The functions J (e) and J'(e) 
are given by the expressions 


J<4V [ 


j xe. 








-CKjp 


A / 
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The second basic integral is of the form 

I- = C 2 /C2 Y 2 ) / M de e‘ (a+2) J(e) 

MAX { e , 2MeV/ 2y } (24) 

for integrations over the Gaussian part of the cross section and a power law 
photon field existing over the energy range e m S G * e M with the form 

n (e) = K e~ a (25) 

The constant in equation (24) has the value 

C 0 = KW*" 1 ^ .S. A -1 (26) 

2 A,i d 

The other two integrals are those over the high energy nonresonant 

part of the cross section above 30 MeV. With the cross section a . = const- 

A,x 

ant over this energy range as given by equation (16 ) , the relevant integrals 
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take the form 


V 


[ j ^ - (z/) j e ^T 


for a thermal photon spectrum n(e) and 


-- K< 5 


MA ‘ 




in the case of a power law photon spectrum. 


Equations (27) and (28) assume a constant cross section for all 
e' above 30 MeV. In reality, photopion production becomes important above 
150 MeV and changes the shape of the cross section. However, because of the 
steeply falling nature of the photon spectrum at high energies, these effects 
are not important for nuclei with Lorentz factors below 10 11 ( energy per 
nucleon below 10^ GeV ). 


Equation (14) was evaluated numerically using the above formulas in 
conjunction with the data given in Tables 1 and 2. Photodisintegration rate 
for nuclei with atomic mass numbers up to 56 were thereby evaluated. (For the 
three nuclei 3 He, 4 He and 9 Be with nonresonant channels having large values of 
A ( A >> 10 MeV ), the 30 MeV cutoff formalism'was not 'employed. )- 

Figures 4 through 6 show samples of the results obtained giving R ft \ 
as a function of y- It is interesting to note how the width of the low energy 
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9 

part of the cross section A affects the rate calculations for y ^ 3 x 10 . 

Another important effect is the difference in the energy dependence between 
the reaction channels dominated by the giant resonance ( i = 1,2) and those 
channels which only contribute to the high energy part of the cross section. 
By combining the rates froija all the reaction channels for a given nucleus , the 
effective nucleon loss rate for the nucleus can be calculated from the relation 


R eff ,A 


(dA/dt) = Z iR . 

i A,i 


(29) 


The quantity R^^fy) is shown in Figure 7 for •' .Fe along with the single em- 

-1 

ission timescale R__ in order to demonstrate the importance of multinucleon 
emission processes in determining the effective nucleon loss rate. 


The different energy loss rates for nuclei can now be compared. Figure 8 

1 56 

the characteristic energy loss timescale t = (dE/dt) -± E for Fe from the ef- 

feet of redshift losses, photodisingtegration and pair production losses. The 

photodisintegration timescale is shown for both the LI R and HIR cases. It can 
56 

be seen that for Fe , the energy loss rates from; pair production and photo- 

20 

disintegration are comparable in the energy range near lCr u eV. 

For photodisintegration, the average fractional energy loss rate 

E“ 1 (dE/dt) = A -1 (dA/dt) (30) 

Jiecause the nucleon emission is isotropic in the rest system of the nucleus. 
The Lorentz factor, which is equivalent to energy per nucleon, is conserved 
fpr photodisintegration as 'opposed to the pair production and photopion prod- 
uction processes which involve the creation of new particles which carry off 


energy . 
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_2 

The energy loss time for pair production t e * AZ as implied by equa- 
tion (4) which increases The energy loss time from photodisintegration 

implied by equation (30) is almost independent of A. Thus, the relative import- 
ance of pair production as an energy loss process vis-&-vis photodisintegration 
decreases with decreaing Z when we consider nuclei lighter than iron. 


5, Monte Carlo Calculations of the Photodisintegration Histories of Nuclei . 


As can be seen from the results of the previous section, the photodis- 
integration rates for all nuclei have typically steep y dependences. Pair 
production can also be quite sensitive to y. Photodisintegration conserves 
y but pair production losses continually and cumulatively act to decrease 
the value of y thereby directly affecting the photodisintegration rate. Thus, 
after a given period of time nuclei with the same initial value of y aquire 
differing new values of y depending on their individual disintegration histor- 
ies A (t) due to the interplay of the differing energy loss processes. The 
only way to take correct account of the effect of changes in y and A' with 
time is to keep track of each individual nucleus as it cascades down in A. 

This is especially important for the low infrared (LIR) case. 


These considerations lead us to use a Monte Carlo technique to compute 
the behavior of nuclei starting at various initial values of y (designated 

by y . At each step, each nucleus followed in the calculation is character- 

0 

ized by the two quantities A^ and (k being the step number). Fofc a given 
time step, the probability of each reaction A^ going to A^-i can be calculated 
using the rates computed previously as functions of y. The step size At is 
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chosen so that the sum of the probabilities for all possible i is smaller than 
l t the difference from unity being the probability for no interactions in a time 
At. Since photodisintegration conserves y, the new value of y is then computed 
for each nucleus using the pair production rate. 

Using the Monte Carlo technique, for each step we randomly determine the 
reaction channel for each nucleus in the system by assigning a random number 
to the nucleus and associating the occurrance of a particular channel with a 
set of such numbers such that the size of the set is proportional to the branch- 
ing ratio for that channel. The Principle of Large Numbers then states that, 
given enough nuclei in the system considered, the frequency of randcm occur- 
rances followed in the calculation will reflect the true probability distribu- 
tion of the reactions as determined by their branching ratios. 

The complete photodisintegration of iron takes, on the average, about 20 
events if one includes multinucleon channels, which indicates that choosing a 
computation run involving a few hundred time steps corresponds to a small in- 
tera±ion probability for each step. The total number of particles in the sys- 
tem must also be such that the number of particles having a value A’ near <A> 
at a given time t is large enough so that a fair representation of the branch- 
ing ratios of different reaction channels can be obtained after they all under- 
go one interaction. 

For each inital value y used in a computation run, 500 particles were 

0 

followed in the photodisintegration cascade. The partial results of two such 

runs are shown in Figure 9 which shows the cumputed histories of two beams of 
56 

Fe nuclei staring out at different initial energies. For comparison, we also 
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56 

indicate in the figure the lifetime for Fe against one nucleon loss for each 

initial energy. These runs were made using the HIR case and it is quite obvious 

that, even in this case, the lifetime for complete photodisintegration of iron 

is ^200 times higher than for the single nucleon loss process. Figure 10 shows 

the computed values of <y{t)> for various values of y for 56 Fe. Also shown for 

0 

comparison is the function y(t) for energy loss of protons from pair production 

for one typical value of y . 

0 

In a few test cases, the same cascade was run several times with a differ- 
ent set of random numbers and different values of At, Typical differences in 
<A> introduced by such changes were of the order of £ 1 percent. This proved that 
the number of particles and time steps used in the computation yielded results 
with adequate precision. 

Using the computed values of the functions <A(t)> and y(t) obtained for dif- 
ferent initial energies, the effect of photodisintegration on a power law initial 
energy distribution can be calculated. This procedure is preferable to the intro- 
duction of a power law distribution in energy of the nuclei used in a Monte Carlo 
run in which case poorer results are obtained because in that case inefficiencies 
would result from there being either too many or too few particles in some energy 
ranges due to the steepness of the initial spectrum. 

For comparison with the results of the Monte Carlo calculation, the quan- 
tity dA/dt for 56 Fe shown in Figure 7 can be used together with the approxima- 
tion that the photodisintegration cross section scales like NZ/A - A/4. Thus, 
for the HIR case for which photodisintegration is significantly more important 
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than pair production 



which implies that A(t) is an exponential function of time with an e-folding 

time of 56(dA/dt^, Figure 11 shows the function <A(t)> obtained from the 

Monte Carlo calculations for various y . for both the HIR and LIR cases. We 

0 

also show for comparison the value of A(t) obtained under the approximation 
given by equation (31). In the HIR case, A is not decreasing as fast as the 
simple scaling approximation (31) would imply because of the diminishing rel- 
ative importance of multinucleon channels for lighter nuclei. For the LIR 
case it is apparent that more than just a change in the photodisintegration 

rate is occurring because of the definite change in the character of the func- 
> 15 

tion A(t) for t n, 5 x 10 s. This change is due to the effect of pair prod- 
uction on y which can be seen to be qualitatively more than a small correction 
and illustrates the rationale for using the Monte Carlo method in order to ob- 
tain results which are truly valid. 3?he shaded areas around the curves shown 
in Figure 11 show the one standard deviation spread in A around <A> determined 
from the Monte Carlo computations from more detailed information as illustrated 
in Figure 9. More <A(t)> curves are shown in Figure 12, By combining these <A(t)> 

results with those on y(t) for various values of y , one can obtain the function 

0 

<A(y^)> where y^ is the final value of y for nuclei starting out as Fe after some 
time t when they have disintegrated to nuclei of mass A. This function for vard.^ 
ous values of t is shown in Figure 13. 


^■The actual dependence of the photodisintegration cross section on A is scmewhat 
higher as can be seen from the values of £ given in Table 1. 
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6. Astrophysical Implications of the Results . 

The results of the Monte Carlo computations discussed in the previous 
section can be ured to discuss the effect of photonuclear emission processer; 
on ultrahigh energy nuclei in intergalactic space by evaluating their effect 
on an initial beam of 56 Fe nuclei with a power law energy spectrum at in- 
jection. One may then ask how this energy spectrum changes after a given 
period of time corresponding to the mean age of the nuclei under various as- 
sumptions as to their origin. The predicted characteristics of the ultrahigh 
energy cosmic rays under these various assumptions may then be compared with 
the data regarding the energy spectrum of ultrahigh energy cosmic rays given 
in terms of total energy per nucleus from air-shower observations. 

We shall consider here two hypotheses of particular astrophysical sig- 
n ^-fi- cance * viz., the universal origin hypothesis and the local supercluster 
or i9i n hypothesis. For the universal origin hypothesis, we can consider that 
cosmic rays have been produced at a roughly uniform rate throughout a time t 

H 

which is related to the Hubble constant H Q and given by (2/3) H* 1 for an 

Einstein-de Sitter universe and H * for a low density universe where H is 

0 0 

taken to be 50 km-s ^*Mpc. The average age of the cosmic rays is then t^/2. 

For the supercluster origin hypothesis, we assume that most of the cosmic 
rays we observe at ultrahigh energy are produced within the local supercluster 
of galaxies which is centered on the populous Virgo cluster at a distance of 
'v 20 Mpc and of which our Galaxy is considered to be a member. Under this hy- 
pothesis, the average age of the ultrahigh energy cosmic rays will be ^xlO 15 s. 
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Of course, an equilibrium spectrum of such cosmic rays produced from a 

uniform distribution of sources will contain a spread in values of A at a 

given value of y. Nevertheless, the functions shown in Figures 11 and 12 

giving A(t) show that a dispersion of a facfcbr of 2 in the age introduced 

a dispersion in A which is not larger than the dispersion introduced by the 

cascade process itself. To simplify the discussion we shall first consider 

the time behavior of a power law energy spectrum of Fe nuclei injected in 

a burst at an initial time t = 0. Starting with such an initial spectrum J(y ) 

0 

and knowing y (t) and A(t), there are then two well defined functions y_(y ) 

f f 0 

and A(y^). This is because the change in y with time is determined by the 
pair production and redshift energy loss rates and the photodisintegration pro- 
cess changes A without changing \ The spectrum in y_at time t is is then ob- 

x 

tained from the inverse function y (y.) using the relation 

0 £ 

I(y^;t) - (Y Q /Y f )d{ ( Y f ) } (32) 

> 

The final energy of the nuclei is, of course, E = y A$-y f )m . E(y) is a non- 

r r p 

monotonic function which reaches a maximum at a critical energy E^, then de- 
creases again to a minimum and starts increasing again with y. This implies 
that nuclei having an energy E = E c + 6E must have started out with an initial 
y much larger than those having an energy E c because of the rapid change of 
A with y from photodisingtegration. Naturally, those nuclei having started out 

with a much higher value of y were produced with a much lower intensity. The 

0 

physical result of this is then that there is a resultant discontinuity or ef- 
fective cutoff in the predicted equilibrium spectrum at an energy E In real- 
ity, the cutoff will not be infinitely sharp because of the dispersion in A as- 
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sociated with the cascade (see, e.g. , Figure 11) and the dispersion in age of 
the cosmic rays due to the spatial distribution of the sources. Nonetheless, 
there will be a sharp drop in the cosmic ray intensity near the value of E c 
associated with the mean age of the cosmic rays. 

Starting with an injection spectrum J(y ) = Y 3 , the resultant spectrum 

0 0 

1 (E;t) was calculated for various values of t associated with the universal 

3 

origin and supercluster origin hypotheses. The function E I(E;t) is plotted 
in Figure 14 for the HIR case and in Figure 15 for the HIR case, given in 
arbitrary units. By comparison, the observational data on the energy spec- 

3 

trum E 1(E) as given by Hillas (1974) are consistent with a straight horizon- 

-3 20 

tal line ( 1(E) E ) up to the highest observed energy of % 2x10 eV. 

The value of E c given as a function of time is shown in Figure 16. 

From this figure, it can be seen that a typical dispersion in cosmic ray 
age of a factor of 2 will produce a change in E c by less than a factor of 

2 so that the predicted cutoff energy as a function of cosmic ray age is 
quite well defined. Thus, the low cutoff energy predicted for the universal 
cosmic ray hypothesis on the assumption that the initfalf cosmic rays are 
intermediate mass nuclei is in clear contradiction with the observations. 

In the case of the supercluster hypothesis however, the predicted cutoff 
energy is too close to the upper limit of the observations, especially con- 
sidering the error involved in determining the energy of a shower-producing 
20 

particle in the 10 eV range, to imply any contradiction. Although the re- 

56 

suits shown in Figures 14 and 15 are given for an initial composition of Fe, 

the values of E for nuclei of lower <A> will be even lower because E “ A for 
c 

a given value of y and the characteristic energy loss rate from photodisinte- 
gration is almost independent of A at a given value of y. 
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The secondary nucleons produced in the photodisintegration reactions, which 
eventually become protons, are emitted isotropically in the rest frame of the 
parent nuclei and thus have, on the average, the same Lorentz factor as the par- 
ent. These secondaries lose energy through pair production and photcmeson prod- 
uction interactions at a rate which is slower than that of the parent nuclei. 

Nevertheless, for a given value of y such that y is a fairly good approxi- 

0 r 

mation for the average y of the secondaries. The spectra of the secondaries ob- 
tained under this approximation are shown in Figures 14 and 15 assuming no cut- 
off in the injection spectrum. For energies above E c , for which the secondaries 

are dominant, the value of y associated with their parent primaries is very 

0 

large and the associated value of A has been 1 for most of the time considered, 
making our approximation a good one. It can be seen that the existence of a sec- 
ondary proton component does not affect our prediction of a sharp cutpff in the 
observed spectrum at E c * 
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7. Conclusions . 

The results of this paper clearly indicate that ultrahigh energy cosmic 

rays therefore cannot be universal in origin regardless of whether they are 

protons or nuclei . (The cutoff energy for protons can be seen from Figure 2 
19 

to occur at ^5x10 eV due to photomeson production losses with a lower en- 

18 

ergy steepening predicted to occur at ^2x10 eV due to pair production, ) . 

Both the supercluster and galactic origin hypotheses, however, are possible 

2 

for protons (Stecker 1968) and higher mass nuclei . 

These conclusions are further supported by reports of a distinct ani- 
' 19 

sotropy in the arrival directions of cosmic rays ,d.th energies above 10 eV 
( Krasilnikov et al. 1974 ) again arguing against a universal origin for 
these cosmic rays since a universal origin implies large-scale isotropy. 

17 19 

When measurements of the mean mass of cosmic rays in the 10 to 10 

eV energy range become available, the results given here for A(y f ) as shown 

in Figure 13 can be used to help determine their age and origin further. 

Should it be found in the future that the cosmic ray spectrum continues to 

20 

values of total energy much in excess of 10 eV, this would rule out their 
being of mixed or intermediate mass composition or of intermediate mass or- 
igin even under the supercluster origin hypothesis. 


2 ! ' ' 

Using new estimates of the infrared radiation field in the inner galaxy, as 

well as estimates of the galactic optical photon field, it can be shown that 
photodisintegration of nuclei travelling typical galactic distances is neg- 
ligible. For a further discussion of the galactic origin hypothesis for ultra- 
high energy cosmic rays, see Stecker (1971) and Syrovatskii (1971). 
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Figure 
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Figure 


CAPTIONS 

1. A graph of the function f (v) used in equation (4) as defined by 
Blumenthal (1970) . 

2. The energy loss time and attenuation length for protons from pair 

production, redshift, and photopion production losses based on the 

calculations of Stecker (1968) and Blumenthal (1970) . Also shown 

56 

is the attenuation length for Fe from pair production losses. 

3. Computed background radiation fields from quasars, Seyfert and 
related galaxies and normal galaxies as discussed in the text. 

Also shown are the 2.7 K microwave and optical radiation fields 
as well as some measured values and upper limits. The dotted 
lines labelled HIR and LIR were taken as alternative models for 
the intergalactic infrared background fields used in the cal- 
culations. 

4. Rates for i-nucleon loss from nuclei of mass number A, labelled 

r" 1 in seconds for an aopropriate sample showing the effect of 
A,i 

different cross section widths A on the shape of R(y) • 

5. more sample inverse rates (see caption for Figure 4). 

6. More sample inverse rates (see caption for Figure 4). 

56 

7. The effective inverse nucleon loss rate for Fe shown along 

with the single emission timescale to demonstrate the importance 
of multinucleon emission processes. 

8. Energy loss time for Fe shown for redshift, pair production, and 

photodisintegration losses in the HIR and LIR cases. 

56 

9. Calculated histories of two beams of Fe nuclei with two different 
initial energies shown for the HIR case. Also shown in both cases 
are ten times the value of the lifetime for Fe against one nucleon 
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Figure 10. 


Figure 11. 


Figure 12. 


Figure 13. 


Figure 14. 


Figure 15. 


Figure 16. 


loss. The gap in the distribution for A between 4 and 9 is caused 

9 

by the reaction y + Be -*• 2a + n. The intensities shown for A = 1 

and 4 are underestimated since they do not include the contribution 

of the secondary products of the photoemission interactions. 

56 

The computed value of y(t) for Fe with given initial values Y q 

shown by the solid lines. The dashed line shows the function y(t) 

9 . 

for protons with y q = 3 x 10 for comparison. 

Sample values of <A(t)> for two values of y • The H1R case is shown 

o 

by the solid lines, the LIR case is shown by the dashed line and 
the scaling approximation (see text) is shown by the dot-dash line. 
More sample values of <A(t)>. Solid lines: HIR case, dashed lines: 
LIR case. 

The function Afy^) calculated for various timescales of astrophysical 
significance. Solid lines: HIR case, dashed lines: LIR case. 

Calculated shape of the ultrahigh energy cosmic ray spectrum after 

_3 

a time t from injection of a y primary Fe spectrum given as 
E^I (E) in arbitrary units for the HIR case and astrpphysically 
significant timescales with the effect of secondary protons included. 
Calculated shape of the ultrahigh energy cosmic ray spectrum in 
the LIR case (see caption for figure 14) . 

Value of the calculated cutoff energy E^ for ultrahigh energy cosmic 
rays starting out as Fe and evaluated after a time t. 
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